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1  Introduction 


This  AASERT  grant  for  the  project  entitled  “Intelligent  Flight  Control  of  Uninhabited 
Aerial  Veiiiclcs”  began  I  June  1997.  This  report  is  the  Final  Report  for  the  period  1  June  1997 
-  30  April  2000.  The  grant  number  of  this  project  is  F49620-97- 1-0406.  The  parent  award  for 
this  grant  w^us  F49620-95-1-0019.  This  grant  ended  30  September  1997.  The  current  parent 
grant  is  F49620-98- 1-0037. 

The  principal  investigator  for  this  grant  is  Professor  Dennis  S.  Bernstein.  This  grant  sup¬ 
ported  the  research  of  Dr.  Tobin  H.  Van  Pelt,  who  completed  the  Ph.D.  degree  in  aerospace 
engineering  at  the  University  of  Michigan  in  April  2000.  Dr.  Van  Pelt  is  a  U.S.  citizen. 

The  goal  of  this  project  was  to  develop,  implement,  and  demonstrate  intelligent  flight  control 
technology.  During  the  previous  two  reporting  periods,  construction  of  the  flight  test  vehicle 
Solus-  was  completed  and  flight  tests  were  conducted  validating  various  aircraft  systems.  The 
paper  [1]  describes  Solus‘S  and  its  research  objectives.  Financial  support  for  the  hardware 
development  of  this  vehicle  was  provided  by  the  University  of  Michigan. 

Theoretical  research  under  this  grant  focused  on  developing  identification  methods  that  were 
suitable  for  on-line  implementation.  The  most  recent  research  reported  in  [2]  concerns  the  de¬ 
velopment  of  a  new  linear  identification  method.  This  research  improves  standard  least  squares 
estimates  by  applying  a  quadratic  constraint  on  the  system  transfer  function  coefficients.  This 
work  is  being  prepared  for  submission  to  IEEE  Transactions  on  Automatic  Control.  Previous 
research  in  [2]  and  [3]  developed  nonlinear  identification  methods.  These  methods  are  least 
squares  biised  and  identify  block-structured  systems  that  contain  static  nonlinearitics  and  dy¬ 
namic  linear  blocks,  and  can  be  implemented  on-line.  The  work  in  [2]  aiul  [3]  has  been  submitted 
to  the  International  Journal  of  Control.  Furthermore,  previous  research  in  [4]  considered  linear 
identification  methods  that  used  sparse  over-parameterizations.  This  research  is  closely  related 
to  adaptive  control  methods  developed  under  the  parent  grant  [5]  and  [6]. 


2  Identification 

Empirical  or  data-based  modeling,  generally  referred  to  as  system  identification,  plays  an 
essential  role  in  control  systems  engineering  as  well  as  many  other  branches  of  science  and 
engineering.  While  advances  in  robust  and  adaptive  control  methods  for  linear  and  nonlinear 
sysUnns  have  helped  to  decrease  reliance  on  plant  models,  the  fact  nmiains  that  models  continue 
to  play  an  important  role  in  controller  synthesis.  Analytical  models  are  essential  for  plant  and 
control  arcliitecture  design  prior  to  plant  fabrication,  but  these  models  contain  assumptions 
that  may  be  very  difficult  or  impossible  to  verify.  Consequently,  siich  models  are  rarely  accu¬ 
rate  enough  to  provide  precise  prediction  of  the  plant  behavior  for  tuning  high  performance 
controllers.  On  the  other  hand,  models  obtained  from  system  identification  methods  incorpo¬ 
rate  the  '‘rear*  system  in  a  more  direct  manner  through  measured  data,  and  thus  reduce  the 
d(q)end(mc('  on  modeling  assumptions. 

Anotlu'r  use  for  systems  identification  is  within  adaptive  control.  Adaptive  control  methods 
r(j(iuir(j  an  accurate  niochJ  of  the  system  in  the  presence  of  changing  disturbances  and  when  the 
syst('m  itself  is  varying  due  to  anomalies  such  as  sub-system  failun^s,  environmental  changes,  or 
system  ch^torioration.  On-line  identification  techniques  provide  models  of  the  plant  in  real-time, 
and  thus  capture  information  about  the  changing  system  to  be  us(^d  by  the  adaptive  controller. 
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One  of  the  first  challenges  of  system  identification  is  the  determination  of  a  model  structure. 
Linear  models  have  been  used  extcmsively  in  control  systems  technology,  but  in  practice,  all  real 
systems  possess  nonlincarities.  Even  when  the  plant  dynamics  can  be  well  represented  by  linear 
models,  it  turns  out  that  additional  nonlinearitics  such  as  saturations  and  dead  bands  can  be 
introduced  by  sensors,  actuators  and  amplifiers,  thus  degrading  the  ciccuracy  of  the  resulting 
linear  model.  Accordingly,  both  liiiear  and  nonlinear  identification  techniques  are  needed. 

An  additional  challenge  of  system  identification  is  the  presence  of  disturbances  that  corrupt 
data  obtained  from  an  identification  experiment.  In  practice,  these  disturbances  could  arise 
from  sensor  inaccuracies,  actuator  granularity  or  ambient  process  noise.  Sometimes  these  dis¬ 
turbances  are  known  and  can  be  removed,  but  in  general  this  is  not  the  case.  Whether  or  not 
these  ambient  disturbances  can  be  removed  during  the  identification  experiment,  it  is  clear  that 
one  of  the  primary  challengers  in  system  identification  theory  is  the  development  of  methods 
that  are  insensitive  to  noise-corru[)t(rd  measurements.  Additionally,  when  considering  nonlinear 
systems  it  is  increasingly  more  difficult  to  distinguish  between  the  presence  of  noise-corrupted 
measurements  and  actual  nonlinear  phenomena. 

One  of  the  most  commonly  used  identification  techniques  is  the  prediction  error  method  with 
a  quadratic  cost  function.  Wlurn  the  model  structure  to  be  identified  is  linear  in  parameters,  this 
technique  leads  to  a  linear  least  squares  problem  which  has  an  analytic  solution  that  is  Ccusily 
implementable.  Our  research  luus  focused  on  identifying  linear  and  nonlinear  model  structur(‘s 
that  can  be  identified  using  a  locust  squares  cost.  Moreover,  we  have  extended  the  standard  least 
squares  theory  by  incorporating  quadratic  constraints.  The  motivation  for  using  le<^ist  squares 
methods  is  that  these  methods  have  closed  form  solutions,  are  easily  implemented  on-line,  and 
continue  to  be  one  of  the  most  commonly  used  mathematical  tools  in  all  of  engineering. 


3  Linear  Identification 


Least  squares  is  a  widely  studied  technique  for  data  analysis  and  system  identification,  and 
it  has  application  throughout  applied  mathematics  and  engineering.  A  critical  issue  associated 
with  least  squares  identification  is  the  effect  of  noise  on  the  parameter  estimates,  A  desirable 
property  of  any  system  identification  procedure  is  that  the  estimates  be  consistent,  that  is, 
the  parameter  estimate  converges  with  probability  one  to  the  true  parameter  as  the  number 
of  data  points  increa.s<\s.  In  the  ca.se  of  least  squares  it  is  well  known  that  the  parameter 
estimates  are  not  generally  consistent.  However,  consistency  is  achieved  when  the  residual 
error  is  uncorrelatcd.  Two  spet'iHc  cases  in  which  this  occurs  are  the  case  of  white  equation 
error  and  the  case  of  a  finite  inqiulse  response  (FIR)  model  structure. 

The  case  of  white  e([uation  en  ror  is  a  highly  unrealistic  situation.  Alternatively,  FIR  models 
can  be  u.scxi  within  least  sf[uan‘s.  In  tliis  case  all  of  the  transfer  function  poles  are  constrained  to 
b(^  zero,  which  is  clearly  a  non-valid  a.ssumptioii  for  infinite  iirqiulse  response  (HR)  systems.  In 
addition,  the  numerator  coc^flicicmts  of  a  FIR  transfer  function  are  given  by  the  impulse  res[)onse, 
or  Markov  parameters  of  tlie  .system.  Hence,  it  is  interesting  to  note  tluit  the  least  scpiarc's 
estimates  of  the  numerator  coefficients  arc  consistent  estimates  of  the  Markov  parameters  of 
the  system  even  if  tlie  actual  syst(nn  is  HR. 

Ther(‘  exist  many  variants  of  the  least  squares  method  that  attempt  to  remedy  the  lack  of 
consistency.  Generalized  locust  squares  simultaneously  builds  a  model  of  the  system  and  tlu^ 
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noise,  and  is  a  special  case  of  tlu*  [)rcdiction  error  method.  Although  generalized  least  squares 
is  consistent  when  using  this  method,  the  legist  squares  cost  is  nonquadratic  in  the  parameters 
and  requires  numerical  optimization  procedures.  In  addition,  consistency  is  guaranteed  only 
for  the  global  ininimizcr  which  may  be  difficult  to  obtain  numerically. 

Another  extension  of  least  squares  identification  is  the  instrumental  variable  method,  which 
uses  constructed  data  sequences,  known  cis  instruments,  that  are  uncorrelated  with  the  noise;. 
It  has  been  shown  that  the  instrumental  variable  method  is  generically  consistent.  However, 
for  finite  data,  the  accuracy  of  the  instrumental  variable  method  may  be  poor  and  in  general, 
the  choice  of  good  instruments  de[)ends  on  the  system  and  the  noise  in  question. 

When  performing  least  sepiares  identification,  a  specific  parameterization  of  the  system 
transfer  function  is  assumed.  In  general,  these  pararneterizations  have  an  inherent  ambiguity 
which  results  from  the  simultaneous  scaling  of  the  numerator  and  denominator  polynomials  of 
the  transfer  function.  It  is  traditional  to  normalize  the  leading  denominator  coefficient  of  this 
transfer  function  to  be  unity  (ao  —  1).  Virtually  every  treatment  of  least  squares  identification 
proceeds  with  the  tacit  assumption  that  this  normalization  entails  no  loss  of  generality  with 
respect  to  least  sc^uares  optimization.  Moreover,  the  transfer  function  is  generally  assumed  to 
be  fully  parameterized.  The  research  in  the  next  two  sections  loosens  these  assumptions  on  the 
system  parameterization  in  order  to  reduce  the  effects  of  noise  and  obtain  consistency. 

3.1  Quadratically  Constrained  Least  Squares 

In  this  research  we  consider  a  more  general  parameter  normalization  for  least  squares  iden¬ 
tification.  Specifically,  instead  of  fixing  uq,  we  enforce  a  quadratic  parameter  constraint  on  all 
of  the  system  parameters.  This  ap[)roach  admits  the  normalization  ao  =  il  as  a  special  case. 

The  least  squares  problem  with  a  quadratic  constraint  was  considered  in  [7].  This  problem 
has  the  form 

min  ||/l:r  - /;||2  subject  to  HCx  —  ri||2  =  a,  (1) 

X 

where  A  €  x  €  ,  C  €  b  G  d  G  W,  and  a  >  0.  In  the  present  research  we 

consider  the  quadratically  constrained  least  squares  (QCLS)  problem 

nnn^^^  A/0  subject  to  9^ NO  =  7,  (2) 

where  M  G  K(2^+2)x(2ri-{-2)  j,.  nonnogativc  definite,  N  G  IR(2^^+2)x(2n+2)  symmetric,  0  G 
and  7  >  0.  In  (2),  0  represents  the  model  parameter  vector 

...  an  bo  •••  6„ 

coiTesponding  to  the  rttli-orcler  transfer  function 

n  -  I  +  •  •  •  + 

ao+a,q-i  +  ---  +  a„q-"’ 

and  9^ MO  is  the  square  of  tin'  residual  model  error.  Note  that  .setting 


is  equivalent  to  enforcing  the  normalization  ao  =  ±1.  Thus,  the  QCLS  problem  specializes  to 
the  standard  problem.  It  can  be  seen  that  (1)  with  6  =  0  and  d  =  0  has  the  same  form  as 
(2)  with  M  =  A  and  N  =  C^C.  However,  for  the  case  of  system  identification,  N  may  be 
indefinite,  which  is  not  allowed  in  the  framework  of  [7]. 

In  this  research  we  show,  rather  surprisingly,  that  the  parameter  normalization  ao  =  1 
entails  a  loss  of  generality  with  respect  to  least  sciuares  identification.  Roughly  speaking,  pa¬ 
rameter  normalization  (however  invoked)  and  least  scpian's  optimization  are  not  commutative 
o[)erations;  in  fact,  each  choice  of  normalization  (or  ([uadratic  constraint  matrix  N)  potentially 
leads  to  a  different  niiniinizcr  of  the  least  squares  cost.  Consequently,  the  normalization  ao  =  1 
actually  subverts  parameter  consistency  except  for  the  white  equation  error  and  FIR  modeling 
Ccises  mentioned  above. 

In  the  present  research,  for  the  case  of  finite  data,  we  show  that  when  a  persistency  condition 
is  satisfied,  there  always  exists  a  quadratic  constraint  matrix  N  such  that  (2)  yields  the  true 
parameters  of  the  system  as  a  solution.  However,  the  appropriate  constraint  matrix  N  depends 
on  the  noise  realization  and  thus  cannot  be  implemented  in  practice.  Furthermore,  when  a 
persistency  condition  is  satisfied,  we  show  that  there  always  exists  a  quadratic  constraint  matrix 
A^,  which  depends  on  the  noise  statistics,  such  that  the  QCLS  estimate  is  consistent.  It  turns 
out  that,  in  the  case  of  white  equation  error  and  FIR  modeling,  the  appropriate  constraint  is 
equivalent  to  fixing  ao-  For  practical  implementation  with  finite  data,  we  develop  an  iterative 
method  for  applying  the  technique  that  does  not  depend  on  knowledge  of  the  noise  statistics. 

Numerical  results  demonstrate  the  improvement  using  QCLS  identification.  A  2nd  order 
example  was  considered  with  colored  noise  resulting  in  a  signal-to-noise  ratio  of  15,  and  a  two- 
tone  sinusoidal  excitation  input.  Figure  1  shows  the  parameter  estimate  error  for  a  typical  run 
using  standard  least  squares,  an  instrumental  variable  method,  the  QCLS  method  with  known 
statistics  and  iterative  QCLS.  Both  iterative  QCLS  and  QCLS  with  known  statistics  exhibit 
less  parameter  error  than  standard  least  sqiuires.  Furthermore,  the  parameter  estimate  error 
of  the  iterative  QCLS  method  is  similar  to  that  of  the  QCLS  method  with  known  statistics. 
Additionally,  the  frequency  response  of  standard  le^ust  squares,  instrumental  variables,  and 
the  iterative  QCLS  method,  are  shown  in  Figure  2.  Again,  the  reduction  in  the  parameter 
estimation  error  of  the  iterative  QCLS  method  yields  the  most  accurate  model  with  respect  to 
the  frequency  response  of  the  .system. 


Method 

I|A0,||2 

Unstable  Models  (%) 

LS 

2.35  ±  .0036 

0 

IV 

1.09  ±2.19 

37 

QCLS 

.281  ±  .020 

2 

N=Q 

.150  ±.0068 

0 

Table  1:  Average  paramct('r  estimation  error. 

Table  1  shows  tlu^  average  AOf  plus  or  minus  oiu'  standard  deviation,  over  100  runs  (X)nsid- 
ered,  as  well  as  the  ruimlx'r  of  unstable  models,  for  ('ac:h  method.  The  results  sliow  that  the 
instrumental  variable  m(‘thofl  is  highly  uncertain,  and  :17%  of  the  identified  models  unsta¬ 
ble.  Moreover,  the  standard  Ic'ast  sc[uares  method  producc^d  models  with  large  bias.  In  contrast, 
iterative  QCLS  produceul  modc'ls  that  had  less  bi^is  than  both  the  standard  least  squares  and 
the  instrumental  variable  rnetliods.  Additionally,  only  2%  of  the  models  were  unstable  when 


using  it(n*ativo  QCLS. 


3.2  /^-Markov  Parameterizations 

Wliilo  a  F\R  iri()(l('l  structure  can  be  used  to  consistently  estimate  Markov  [)aramcters  of  a  lin¬ 
ear  system,  it  turns  out  that  convergence  can  be  slow  and  more  rapidly  converging  algorithms 
are  of  interest.  One  idea  for  improving  the  speed  of  convergence  is  suggested  by  the  fact  that 
the  leading  nmrKuator  coefficient  of  a  fully  parameterized  transfer  function  is  also  a  Markov 
parameter  of  the  system.  Therefore  it  is  natural  to  seek  model  parameterizations  whose  nu¬ 
merator  polynomials  include  an  arbitrary  number  of  Markov  parameters  without  incorrectly 
constraining  the  pole  locations  of  the  model  structure.  The  objective  in  using  such  models  is  to 
obtain  consist(mt  (‘stimates  of  the  system  Markov  parameters  with  a  faster  rate  of  convergence. 

The  significance  of  this  result  arises  from  the  fact  that  the  estimates  of  the  Markov  param¬ 
eters  can  b('  extracted  and  used  to  construct  a  model  of  the  system  using  realization  theory. 
Since  this  [)rocedur(^  can  be  carried  out  whether  the  system  is  FIR  or  IIR,  it  can  thus  be  seen 
that  an  ‘'incorrect"  model  structure  can  effectively  be  used  to  construct  a  suitable  identified 
linear  model. 

A  suitable  model  structure  that  satisfies  the  above  criterion  is  given  by  the  class  of  step- 
ahead  predictors,  or  equivalently,  ARMARKOV  (ARMA  -h  Markov)  models,  which  are  non- 
minimal  parameterizations,  some  of  whose  numerator  polynomial  coefficients  coincide  with 
Markov  {)arameters.  These  models  have  been  used  extensively  for  adaptive  control  and  identi¬ 
fication  [5-6]. 

Our  prc'viously  reported  research  systematically  developed  and  explored  the  properties  of 
transfer  function  parameterizations  that  explicitly  contain  Markov  parameters  in  the  numerator 
polynomial.  An  interesting  feature  of  these  models  is  the  relation  between  the  number  of 
Markov  parameters  appearing  in  the  numerator  and  the  sparse  structure  of  the  denominator 
polynomial.  VVe  characterized  this  relationship  by  defining  /i-Markov  parameterizations,  which 
arc  non-minimal  parameterizations  of  a  transfer.  The  [)rincipal  result  of  this  research  guarantees 
that  when  white  noise  excitation  is  used,  the  least  sciuares  estimates  of  the  Markov  parameters 
are  consistent  under  a  very  general  class  of  noise  models,  and  regardless  of  model  order  choice. 
These  Markov  parameters  can  then  be  used  with  realization  theory  to  obtain  consistent  model 
estimates.  The  significance  of  this  result  for  obtaining  consistent  model  estimates  is  that  precise 
knowledge  of  tlie  noise  and  the  system  order  is  no  longer  needed. 


4  Nonlinear  Identification 

As  [ireviously  reported,  we  have  considered  nordinear  identification  using  a  model  represen¬ 
tation  involving  linear  dynamic  blocks  (£)  and  static  nonliiu'ar  blocks  {Af).  These  elements 
can  be  interconnected  in  various  configurations,  thus  r('[)resenting  a  large  class  of  nonlinear  sys¬ 
tems.  The  most  fr(H[uently  used  structures  are  tlie  HamnunstcMn  model  {Af  ->  £),  the  Wiener 
model  (C  Af).  and  the  nonlinear  feedback  rnocUfi  (C  Af  £),  while  other  configurations 
are  possible.  Much  of  the  previous  work  in  nonlinear  syst(nns  identification  using  this  model 
structun'  has  assunuxl  that  one  or  more  of  these  sub.systems  are  known  a  priori.  Otlier  work 
has  relaxed  this  ^ussumption  while  using  iterative  procedure's. 

In  our  research  we  develop  a  novel  technique  for  nonlinear  system  identification  using  the 
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Harniuonstein,  nonlinear  feedback,  and  Hairiinersticn/nonlinear  feedback  nioch^l  st,ructure.  This 
technicinc  siinnltaneously  identifies  the  C  and  M  blocks  without  a  priori  inforniation  of  the 
system.  The  key  to  tins  research  is  tlie  development  of  a  parameterization  of  the  nonlinear 
static  maps  that  k^ad  to  a  nonlinear  least  squares  cost  that  can  be  reduccHl  by  means  of  an 
overboundin^  technic[!ie.  This  overbound  entails  a  suboptirnal,  but  computationally  tractable, 
approximation  to  tlu'  original  nonlinear  least  squares  cost.  In  particular,  this  approach  replaces 
the  nonlinear  least  s([uar('s  optimization  problem  witli  a  pair  of  standard  computational  proce¬ 
dures,  nanu'ly,  liiu'ar  k^^ust  squares  optimization  and  fixed  rank  approximation  in  the.*  Frobeuius 
norm. 

This  approach  is  l)ased  on  a  parameterization  of  the  nonlinear  static  maps  J\f  that  is  linear 
in  parameters.  Although  there  are  many  such  parameterizations,  piecewise  linear  functions  are 
a  prime  candidates  since  they  can  approximate  a  large  class  of  nonlinearities.  These  approxi¬ 
mations  have  been  widely  studied  in  the  circuit  theory  literature  as  well  as  in  nonlinear  control 
theory. 

One  of  tlu'  chalksnges  of  this  approach  is  in  guaranteeing  the  invertibility  of  tins  regression 
matrix  that  ari.ses  in  the  linear  least  squares  optimization  problem.  It  is  well  known  that  the 
singularity  of  the  r(*gression  matrix  for  problems  involving  a  linear  model  structure  is  linked  to 
persistency  of  ('xcitation  conditions  on  the  input.  A  key  feature  of  our  method  is  the  develop¬ 
ment  of  a  '‘point-slope"  re[)resentation  for  the  parameterization  of  the  approximating  piecewise 
linear  function.  This  parameterization  guarantees  the  invertibility  of  the  regression  matrix.  VVe 
have  used  this  technique  to  develop  nonlinear  identification  methods  for  Hammerstcin,  non¬ 
linear  feedback,  and  llammcrstein/nonlinear  feedback  model  structures  and  w('  are  currently 
working  on  VVieruM'  and  other  more  general  configurations.  Figure  3  shows  an  (example  of  an 
identified  nonliiu'ar  map  corresponding  to  an  input  saturation.  The  u.se  of  this  repiTscntation 
for  the  nonliiu'ar  static  map  M  is  essential  since  using  other  representations  such  as  standard 
interpolation  n'siilt  in  singular  regression  matrices. 


5  Status  of  Flight  Experiments 

The  goal  of  tlu'  UAV  program  at  the  University  of  Michigan  is  to  develop  and  demonstrate 
intelligent  flight  control  techniques  that  are  relevant  to  UAV  applications  and  advance  the  state 
of  the  art  in  tlu^sc^  areas.  Specifically,  the  research  and  experimental  focus  of  our  project  is  in 
the  fabrication  of  a  tc'stlxxl  that  is  capable  of  performing  on-line  identification,  fault  detection 
and  reconfiguration  of  flight  controls.  The  expected  result  is  a  reduction  in  tlie  risk  and  length 
of  time  associated  with  transitioning  innovative  research  to  flight  control  a[)plications. 

Our  first  aircraft  was  the  UAV  Solus.  Solus  was  a  |  scale  Citabria  with  an  emluKlckxl  control 
system  running  a  r(‘al-time  operating  xSystern  capable  of  communicating  with  a  remote  ground 
station  during  opc'ration.  During  flight  tests  of  Solus,  many  of  our  avdonics  and  telemetry 
systems  were  t(*st(‘d.  During  these  proof-of-concept  flights  it  was  discoverf'd  that  tlu^  Citabria 
airframe  was  not  id(‘al  for  our  application  due  to  weight  and  volume  (*onstraints  and  it  was 
determined  that  a  lu'w  airframe  was  needed. 

Design  and  (‘oust  i  iiction  of  a  new  aircraft  Ix^gan  in  January  of  1998  and  Soliis*^  was  completed 
in  December  of  1998.  JJiis  new  aircraft  us(^s  an  upgrackxl  embedded  flight  control  .system  derived 
from  Solus,  Soln.s-  ntilizc's  a  dual-boom,  pusher  configuration  with  a  5  h.[).  two  stroke  engine. 
The  aircraft  weighs  55  lbs.  and  has  a  131  in.  wings{)an.  A  photo  showing  the  lU’w  aircraft  is 
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Figure  3:  Estimate  of  a  saturation  nonlinearity  in  a  Hainmerstein  system. 


show  in  Figure  5. 

Solus^  carries  an  extensive  package  of  sensors.  It  has  an  inertial  measurenuMit  unit  (IMU), 
an  air  data  system  (ADS),  and  a  Trimble  GPS  receiver  with  an  Accupoint  Dif[er(mtial  GPS 
correction  unit.  The  design,  construction,  and  calibration  of  the  IMU  and  ADS  were  performed 
at  the  University  of  Michigan  over  the  past  year.  Additionally,  there  is  an  engiiu^  tacliomoter  to 
measure  engine  speed  and  potentiometers  mounted  on  each  control  surface  to  measure  control 
surface  deflection. 

The  air  data  system  measures  the  angle  of  attack,  the  angle  of  side-slip,  and  t\\v,  dynamic 
pressure  of  th('  vehicle.  The  angles  are  sensed  with  vanes  connected  to  potentiouu'ters  that 
align  themselves  with  tlie  local  airflow,  and  the  dynamic  pressure  is  sensed  with  a  F^itot  probe 
that  is  connected  to  a  pressure  transducer.  The  resulting  system  is  capable  of  providing  the 
airspeed  of  the  veliick^  as  well  cus  the  wind  direction  during  flight. 

The  iiKutial  m(‘asur(*m(uit  package  consists  of  6  Analog  Devices  solid  states  accehuonieters, 
3  British  Acuospace  solid  state  rate  gyroscopes,  and  a  Honeywell  3  axis  solid  state’  Huxgate 
magnetometer.  Gombiiu’d  with  the  other  sensor  packages  and  appropriate  filtering  the’  complete 
state  of  the  aircraft  can  be  measured.  The  calibration  was  carried  out  at  tlu’  University  of 
Michigan  using  an  Ideal  Aeuosmith  rate  table  and  Singer  Scorsby  table. 

The  UAV  can  l)e  controlUxl  using  two  processors  on  the  aircraft  that  communicati’  through 
dual  ported  R  AM  and  are  in  contact  with  the  instrumentation  through  a  PC  KM  bus.  Addition¬ 
ally,  the  aircraft  is  in  contact  with  a  ground  station  that  consists  of  a  single  120  MHz  Pemtium 
laptop.  This  ground  statioii  provides  the  user  interface,  data  storage,  and  flight  [)lanning  ca¬ 
pabilities.  The  flight  c.om[)utcr  and  ground  station  operate  under  the  QNX  real-tiuu^  operating 
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